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Source modelling 



Forward problem 

bj = J g L^(r) ■ y(r)dV 

V t = lV t (r)V(r)dV 

Inverse problem 

bj(t),V k (t), L"(r),I*(r) -> J p est (r,t) 
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Bayesian approach: basics 


fpoa ( x I b) = f prior (x)/(b I x) 

x the model parameters 

b the measured data 

fprior ( x ) a priori knowledge about the sources 

f post (x I b) probability distribution of the sources after a measurement 

/(b I x) probability distribution of the data for given sources 


MSH 8/2000 



Bayesian probability densities 


• /priori*) incorporatesthea prioriassumptionsabout 
the sources 

• /(blx) tells what kind of data we expect when 
the source fixed, i.e., itcontainsthe solution ofthe 
forward problem forourmodel 

• f post ( x I b) is the 'solution 1 ofthe inverse problem, 
i.e., what is our knowledge about the source 
distribution after the measurement has been 
made 
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Howto characterize the result? 


•Which source para meters maximize f post (x\b)7 
'most probable source 1 , least-squares fit 

• How to choose a set of source parameters 
to satisfy: 

| fpost (xlb)dx = p/100 

xeS 

(p-percent confidence interval, not necessarily 
simply connected: 'local minima') 

•Characterizing a the dependence on a multidi¬ 
mensional source may be difficult 



Why is this helplul? 


•The underlying modelling assumptions 
are explicitly specified 

•The desired quantification of the 
a posteriori distribution is clearly stated 

•The technical problem of obtaining the result is 
separated from the specification of the model 
assumptions 

• It is in principle easy to modify the modelling as¬ 
sumptions (calculation maybe difficult!) 
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Current dipole model 


The neural currents on 
a few-cm 2 patch of cortex 
are approximated 
by a currentdipole 


Time-varying dipole model: 
•dipole positions a re fixed 
• amplitudes vary with time 
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# Model for measured data 1 


Modeldata +additive zero-mean Gaussian noise, 
known comelationsacnossmeasurement points 


Noise 


B = GQ + N^g>,,e,) q ; + N 


Unit-dip oleg a 


Dipole positions 

Dipole orientations 

Dipole waveforms 
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w Model for measured data 3 


No time correlations: 

Einy^Sjfi 

Mean and variance of B: 

E{bj} = (GQ)j 

E{(b j -E{b j })(b k -E{b k }f} = d jk C 

Ateach time pointthe data are jointly Gaussian with 
the above averagesand correlation matrix 



Non-dipolarfields 


Any field patternsg k can be included: 

• Diffuse (distributed) sources 

•Fields of stereotypic environmental noise sources: 
determined from empty room measurements using 
an SVD of the data matrix 

• 'Noise' sources in the subject: MOG, MCG 


lime dependencies of the distributed 
and noise sources can be obtained 
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# Probability densities? 


• Data distribution fora given model: 

^ k 

•A priori distribution: f prior (x) = p pr in the brain, 

0 otherwise 

•A posteriori distribution (0 outside the brain): 

fpo St (x\B) = Ppo J(B\x) 

We have now 'solved'the inverse problem under 
these assumptions! 


/(BI x) = p 0 exp 
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# Most probable source? 


Minimize 

'£(b t -(GQ) 1 ) I 'C- 1 (b t -(GQ) 1 )' 

V k J 

over 

r k AQ 

This is a (weighted) least-squares fit! 

Amplitudes: easy forgiven positions 
and orientations 

Positions: difficult, nonlinear minimization 
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# Extensions a nd modific ations 


• Exponential distribution of noise: robust estimation, 
guard against outliers in the data (LI norm in data 
space) 

• Source correlation matrix: 

a priori distribution forsource amplitudes must be 
modified 
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Dipole modelling in practise 


Visual inspection: 

• Plana rg radio meter signal distribution 

• Normal field component distributions 

Fitting: 

• Single timepoint fits for isolated sources 
•Channel selections 

• Multidipole search with optionally fixed dipoles 
•Alternative: global optimization 



Practical measurement system 



MSH 8/2000 


Bifocal sensors 
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Median newe SEFsignals 
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# Adopted fitting strategy 


1. Fit one dipole at the 20-ms peak. 

2. Fit next dipole at the 35-ms peak. 

3. Using a selection of channels above one of the 
temporal areas, fit the Sll activity on one side. 

4. Use the 20-msand 35-msdipolesin the multidipole 
fitting with fixed positionsand orientations. 

5. Fit the Sll dipoles using the dipole on one 
hemisphere and it's mirror image asinitial 
guesses. 


MSH 8/2000 



Source waveforms and locations 
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w Conclusions 


• The Bayesian approach providesa framework to state 
the modelling assumptions clearly 

• Assumptions can be easily modified 

• Bayesian approach can be employed with all source 
models 

• Numerical difficulties a re clearly separated 

• Heuristic analysisstategiesane still often 
needed with the current dipole model 



